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We evaluate the frequencies of scissors modes for density and concentration fluctuations in a 
vapour of fermionic atoms placed in two hyperfine levels inside a spherical harmonic trap. Both the 
superfluid and the normal state are considered, with inclusion of the interactions at the random- 
phase level. Two main results are obtained: (i) the transition to the superfluid state is signalled by 
the disappearance of soft transverse modes of the normal fluid in the coUisionless regime and (ii) the 
eigenfrequency of the density fluctuations in the superfluid coincides with that of the normal fluid 
in the coUisional regime. The latter property is related to the opening of the gap in the single-pair 
spectrum. 
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I. INTRODUCTION 



The possibility of obtaining a novel superfluid system by pairing of ultra-cold fermionic alkali atoms seems to 
become closer to experimental realization [ l],^ . Various mechanisms of pairing have already been proposed [^-^ , the 
simplest one being an s-wave pairing between atoms belonging to two different hyperfine levels of a magnetic trap. 

Revealing the superfluid transition in a Fermi gas is expected to be more difficult than for bosonic condensation, 
since it cannot be inferred from the observation of the density profile of the cloud . One has therefore to look for 
the transition in dynamical or kinematical properties of the fluid. We suggest in this work a possible method for the 
detection of the superfluid character of the transition. The same test has already been successfully employed in the 
case of a Bose- Einstein condensate [^,^ . 

The property of superfluidity in a vapour is related by definition to the character of its dynamical response to a 
long- wavelength transverse probe which will excite only the non-supcrfluid component. In an inhomogeneous finite 
system, such as an atomic vapour in magnetic or optical confinement or an atomic nucleus, a test of the superfluidity 
can be obtained from the excitation of a small-angle oscillation in a plane where the confinement is slightly anisotropic. 
For a Bose-condensed gas this scissors mode is predicted to show only one frequency component in the superfluid 
state, while it has two frequency components in the normal fluid due to the additional contribution of transverse 
excitations. The study of the scissors mode in Bose-Einstein condensates of alkali vapours has been suggested by 
Guery-Odelin and Stringari Q and an experimental realization of their ideas has already been given by Marago et al. 

0- 

The low-energy collective excitations of a trapped Fermi gas in the superfluid state have been investigated by 
Baranov and Petrov |^ in the dilute (non-interacting) limit. Their theory is formulated in terms of the fluctuations 
in the phase of the order parameter, but they show that these modes also manifest themselves as density fluctuations 
in the sample. Their equation of motion for the density fluctuations of a superfluid Fermi gas coincides with that 
obtained by Amoruso et al. [lO| ] for a non-interacting normal Fermi fluid in the coUisional regime, within a local-density 
approximation (LDA). 

In this paper we derive in Sect. |l| an equation of motion for the density fluctuations of a superfluid Fermi gas in 
an improved LDA, which includes also the effects of the interactions within a random-phase approximation (RPA) 
as already proposed in early work by Anderson jl^. We make use of this equation in Sect. [II to give theoretical 



predictions for the excitation frequency of the scissors mode in the superfluid Fermi gas. This result is compared in 



Sect. IV with the excitation frequencies of the scissors modes of a normal Fermi fluid in the degenerate regime and in 
the presence of interactions. Finally, Sect. ^ presents a summary and some concluding remarks. 



II. COHERENT DYNAMICS OF A SUPERFLUID FERMI GAS 



The low-energy excitation spectrum of a neutral Fermi gas with weak attractive forces at very low temperatures is 
characterized in the superfluid phase by a longitudinal collective mode with a linear dispersion relation at the velocity 
of ordinary (first) sound. This mode is not allowed to propagate in the normal phase as it is embedded in the particle- 
hole continuum of excitations (just as it is the case for the zero-sound mode in a fluid with attractive interactions). 



1 



However, in the superfluid phase the opening of the gap A shifts the energy threshold for the quasiparticle-quasihole 
continuum to 2A, leaving a window in energy where first sound is stable. This mode is known as the Bogolubov- 



Anderson sound |11 1^ and can be viewed as a coherent oscillation of the order parameter for the condensed phase; 
in fact, its presence is required by the Goldstone theorem to ensure gauge invariance. The simplest theory which 
includes the Bogolubov- Anderson mode in the density fluctuation spectrum is the random-phase approximation for a 
superconductor |pT| , |l3| , ^ . 

In a harmonically trapped superfluid the Bogolubov- Anderson sound becomes quantized due to the confinement. In 
the dilute limit the equation of motion for the density fluctuations has been obtained by Baranov and Petrov |^ using 
a semiclassical solution of the Bogolubov - de Gennes equations. As already noted in Sect. |, the resulting spectrum of 
excitations coincides with that obtained ||l0[| for an harmonically trapped normal Fermi gas in the coUisional regime, 
when a local-density approximation for the kinetic stress tensor is employed. 

The coincidence of the two spectra can be explained - within the LDA - by making use of the result for the 
homogeneous fluid that the compressibility is unaffected by the phase transition [ p5| : a superfluid and a normal fluid 
at the same density and interaction strength show the same value for the compressibility. Indeed this property is 
not restricted to the non-interacting Fermi gas: the coincidence of the sound velocity in a superfluid and in a normal 
fluid in the coUisional regime has also been demonstrated by Larkin and Migdal [l^ for generic interaction strengths 
within the theory of superfluid Fermi liquids. For the speciflc case of a weakly interacting Fermi gas the result for the 
sound velocity c in the normal fluid in the coUisional regime, given by = (v'j/3){l + 2fc/a|j^/7r) coincides with the 
velocity of the Bogolubov- Anderson mode, as obtained from the RPA |0. Here, Vf and kf are the Fermi velocity 
and wavenumber, while an is the s-wave scattering length between fermions in different hyperfine states. 

Another characteristic of the spectrum of a homogeneous Fermi superfluid is that all collective modes other than the 
density fluctuations, such as spin-density and transverse current-density modes, are suppressed in the zero temperature 
limit, since they are related to the motion of the normal component of the fluid [ |l8yi9[| . 

These observations allow us to conclude that the equations of motion for the total density fluctuations in a symmetric 
two-component normal Fermi gas in the coUisional regime ||2^ describe as well the collective excitations of a weakly 
interacting superfluid at T = within a local-density approximation. Specifically, we take as the basic equations for 
the dynamics of a superfiuid Fermi gas in the linear regime the following: 

dtn{r, i) = -V ■ {neq{r)v{r, t)) (1) 

and 



9tv(r,t) = -V 



iA(ne,(r)/2)-i/3,,(j., t) + ^gn{r,t) 



(2) 



Here A = /i^(67r^)^/'^/2m, g = AiTh^a-\i/m, and n{r,t) and v(r, i) are the total density fluctuation and velocity field 
for the gas, Ueqij) being the equilibrium density profile. Interactions have been included in the theory at mean field 
(Hartree) level and exchange is not allowed since s-wave interactions are active only between fermions with different 
spin polarization. In the homogeneous fluid this approximation on the interactions leads to the RPA expression for 
the sound velocity [0. The results of Baranov and Petrov |^ are recovered by setting g ~ in Eq. 

It is interesting to notice that Eq. (||) for the velocity field is irrotational, in analogy with the corresponding equation 
for a Bose-condensed gas |2^ . In this case it contains explicitly the effect of Fermi statistics in the expression for the 
linearized local chemical potential. Equations (|l|) and (H) apply only for describing low-energy excitations at very low 
temperature, since they do not include the coupling with thermal quasiparticle excitations and the related damping. 

As an application of Eqs. we study the possibility of exciting the scissors mode in a superfluid Fermi gas. 



III. SCISSORS MODE IN A SUPERFLUID FERMI GAS 



We consider an atomic cloud confined inside a slightly anisotropic trap in the xy plane, described by the confining 
potential 

Vext{r) = + s)x^ + imw^(l - e)y^ + ^mwfz^ . (3) 

Within this geometry the scissors mode can be excited by a sudden rotation of the trap through a small angle in the 
xy plane. In the experiment by Marago et al. the anisotropy has been obtained by adding a small component in 
the z direction to the magnetic field of the TOP trap and the scissors mode has been excited by a sudden change of 
the sign of such field. 
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For the study of the scissors mode of a superfluid Fermi gas we employ the technique of moments developed in |2^] 
(see also We define the dynamical average of a variable x(r, v) on the total density profile n{r,t) as 

(x) = J dVx(r,v)n(r,i) . (4) 

For the scissors mode the relevant variable is x = xy, describing the quadrupolar oscillation with z-component of the 
angular momentum rriz = ±2. By making use of the equations of motion for the superfluid we obtain 



from the continuity equation and 



dt{xy) = {xvy + yvx) (5) 



dt{xvy + yv^) = ~2Lol{xy) (6) 



from the linearized Euler equation. A physical picture of the scissors mode can be obtained by observing that the 
moment {xy) is related to the angle 9 of oscillation of the cloud in the xy plane, \i 6 <^1. For larger angles instead 
the variable {xy) describes the usual quadrupolar excitation of the cloud. 

In deriving Eq. we have made use of the expression for the equilibrium density profile in the local- 
density ("Thomas-Fermi") approximation ]23| , p^ , which for a symmetric two-component fluid is the solution of 
A[ne<j(r)/2]^/'^ = /i — Vext{i^) + 3"-eg(r)/2. We are here exploiting the fact that even in the presence of weak in- 
teractions the equilibrium density profile of the superfluid is well approximated by the corresponding profile of a 
normal fiuid at the same number of particles. This property is illustrated in Figure |l] and is a consequence of the 
interplay between the external confinement and interactions. While interactions shift only slightly the chemical poten- 
tial of the homogeneous superfluid relative to that of the normal fluid, in the presence of confinement they also enter 
in the density profiles through the position-dependent Hartree term. This second effect is dominant in determining 
the profiles, as shown in Figure 0. As for confined Bose condensates, the fluid is dilute but nevertheless interactions 
significantly modify the density profile. 

The solution of the coupled equations (||) and (^) yields for the frequency of the mode the result 

= 2ujI . (7) 

This value coincides with that for the (n = Z = 2) total-density surface mode of a two-component normal Fermi gas 
in the coUisional regime . The same property holds for a Bose-condensed cloud Q . 

In fact, the scissors mode is not affected by the interactions: its frequency coincides with that of the (n = / = 2) 
surface mode in a single-component Fermi gas 10 . This is indeed a property of all surface modes of total density 



fiuctuations |gO| . For these modes V ■ v = and the equation of motion for the velocity field reduces to 

ma2v(r, t) = V [v(r, t) ■ V[A(ne,(r)/2)2/3 + gn,,(r)/2]] . (8) 
By employing the Thomas-Fermi form of the equilibrium profile we have 

mdfv{v, t) = -V [v(r, t) ■ V^xt (r)] • (9) 

This equation shows that only the external confinement determines the frequency of these modes. Equation (^ is 
independent of the statistics: it describes as well the surface modes of a classical gas in the hydrodynamic regime and 
of a Bose-Einstein condensate at T = [Esll • 



IV. SCISSORS MODES FOR A NORMAL FERMI GAS 

We shall now contrast the result (0) for the superfiuid with the excitation frequencies of the scissors modes in a 
two-component Fermi fluid in the normal state. At a temperature higher than the superfluid transition temperature 
Tsup the gas is described by the Vlasov-Landau equation for the Wigner distribution functions /p (r, t) of each of the 
components: 

(St + ^ • Vr - Vr(K.t(r) +gn^(r,i)) • Vp) /^(r,i) = I,M[f^{r,t)] . (10) 

The density profile of each component is defined in terms of the Wigner function as na-(r,t) — J dPpfp{r,t)/{2Trh)^ 
and Icoii is the collision integral. 
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The average of the variable x(r, v) on each component cr of the gas is defined as 

ix). = J d\ J d^px{v.^)f;{v,t) . (11) 

The equation of motion for {xy)^ is now coupled to both the longitudinal and the transverse excitations and to the 
kinetic tensor fluctuations, yielding the following set of equations: 

dt{xy)a {xvy + yvx)a , (12) 

dt{xvy - yv,), = ~2ewl{C{xy), + (1 - C){xy)s) , (13) 

dt{xvy + yv^)„ = 2{W^y)Jm - 2ujl{C{xy)„ + (1 - C){xy)s) (14) 

and 

dt{^xv)a = -mCuj^ [{xVy + yVx)a- + e{xVy - yVx)„] - {Xlxy)a/T . (15) 

We have defined the current density Ua-Va ~ J (Pppfp{r,t)/m{2nti)^ and the xy component of the kinetic stress 
tensor 11!^^^ = / (fppxPyfp{r,t)/m{2TTh)^, and made use of the Thomas-Fermi Ansatz An'^^{r)^^^ = C{Ef — Vext(r)) 
for the equilibrium density profile in a symmetric system, as already employed in The constant C — (Ep/ Ep)'^ , 
where Ep is the true Fermi energy and Ep is that of the non-interacting Fermi gas, measures the strength of the 
interactions in its deviations from unity. The approximations that we have made require |C — 1| <C 1 and T <g; Tp; 
the latter condition is well compatible with the condition T > T^^p in a weakly interacting fluid. Notice that the 
motions of the two components of the gas are coupled for C ^ 1 . 

Collisions have been included in Eq. (|l5|) within a single-relaxation-time approximation by setting 

{nxyIcoll)a = -{nxy)a/T . (16) 

This approximation subsumes damping by scattering against impurities and by inter-species scattering (see ]20| ] for a 
full discussion). 

In the collisionless regime {lot » 1), equations (|T^p^ can be combined to yield the following coupled differential 
equations for {xy)^: 

{df + ACuldl + Ae^wtC^){xy), + 2loI{1 - C){-df + 2e^ul){xy)^ = . (17) 

Solution by diagonalization yields two scissors modes associated with total density fluctuations and two further modes 
associated with concentration fluctuations, at frequencies given by 

„2_. 2 r 2Ct(1-C) , 
w - ^'^(iye^c{C±{l-C))/{2CT{l-C)) 

The appearance of soft modes with a frequenc y p roportional to e is a peculiarity of a non-superfluid system. 
In the collisional regime (wt ^ 1), Eqs. ( p^ - |l^ ) instead yield 

(a^ + 2Culdt){xy)„ + 2(1 - C)uoldt{xy)^ = . (19) 

This has two solutions for non-zero frequencies of the scissor modes at 

uj^ ^ 2ujI{C ± {1 - C)) . (20) 

These agree with the (n = / = 2) surface modes already obtained by a different approach in [^. As expected, the 
frequency of the total-density mode reproduces the result in Eq. (|^). 



V. SUMMARY AND CONCLUDING REMARKS 



In this work we have extended the equations of motion for density fluctuations in an inhomogeneous Fermi gas in 
the superfluid state at T = to include interactions at a mean-field (RPA) level within a local-density approximation. 
In the appropriate limits our approach reproduces Anderson result pT| , p^ for the speed of the Bogolubov- Anderson 
phonon in the homogeneous superfluid and those by Baranov and Petrovp] for a confined non-interacting superfluid. 
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As an application, we have computed the frequency of the scissors mode for a fermionic superfluid gas. By comparing 
the result with those for a normal Fermi gas in the coUisionless regime, we have shown that the measurement of the 
scissors modes could be used as a signal for the superfluid transition, since well below the transition temperature the 
soft transverse modes are suppressed. However, these modes are already non-propagating in a normal Fermi fluid in 
the collisional regime. 

We have also found that the excitation frequency of the scissors mode in a superfluid Fermi gas coincides with 
that of a normal Fermi gas in the collisional regime. This result is understood in the fermionic case as related to the 
opening of the gap which stabilizes the first-sound mode. As is the case for all surface modes, interactions do not 
shift the frequency of the scissors mode in the superfluid. 

While we have worked here within a local-density approximation, further work should be addressed to treat in a 
fully quantal way the dynamical properties of the inhomogeneous fluid of present interest. 
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FIG. 1. Equilibrium density profile at T = in spherical harmonic confinement (with trapping fi-equency uJho = in x 520 
s~^) for a superfluid Fermi gas (solid line) and for a normal Fermi gas (dot-dashed line) in the symmetric system at the same 
number of particles A'^ — 6.6 x 10'*. The two curves arc the result of a local-density approximation using the chemical potential 
of the homogeneous fluid, as evaluated in the first case from the solution of the BCS equations, and in the second from the 
non-interacting value iJ,hom[n] = h^{3n^n)^^^ /2m. The equilibrium density profile of a non-interacting normal Fermi gas at the 
same N is also shown (dashed line). 
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